THE CHEBOTAREV INVARIANT OF A FINITE GROUP 



EMMANUEL KOWALSKI AND DAVID ZYWINA 

Abstract. We consider invariants of a finite group G related to the number of random 
(independent, uniformly distributed) conjugacy classes which are required to generate it. 
These invariants are intuitively related to problems of Galois theory. We find group-theoretic 
expressions for them and investigate their values both theoretically and numerically. 



1. Introduction 

A well-known method to compute the Galois group H of a number field (e.g., of the 
splitting field of a polynomial P G Z[T] with integral coefficients) can be described roughly 
as follows: (1) find a group G which contains H, e.g., because of symmetry considerations 
(2) this group G is our "guess" for H, and we try to prove that H = G by reducing modulo 
successive primes, using the fact that the Frobenius automorphisms give conjugacy classes 
in the Galois group H, and hence conjugacy classes in G. 

If the guess in (1) was right, and if one is patient enough in (2) that the conjugacy classes 
observed are only compatible with the Galois group being our candidate G, then we have 
succeeded. 

This method is particularly simple when G is "guessed" to be the symmetric group acting 
on the roots of the polynomial: it can lead quickly to examples of equations with this Galois 
group. In general, however this is not the most efficient algorithm (if only because the first 
step (1) is hard to formalize!), and thus computer algebra systems use other techniques. Still, 
this method is well-suited for certain theoretical investigations, for instance for probabilistic 
Galois theory (see, e.g., [G]), and it can be surprisingly efficient even for fairly complicated 
groups (see for instance our joint works |JKZ] and [JKZ2j with F. Jouve, for a case where 
the Weyl group of the exceptional Lie group E s is the Galois group involved, and for for 
further developments along these lines). 

In view of this, it is somewhat surprising that no general study of the efficiency of the 
underlying algorithm seems to have been performed. Among the very few references we know 
is a paper of Dixon |Dlj . who considers informally the case of the symmetric groups & n and 
mentions some earlier work of McKay. On the other hand, there has been a fair amount of 
interest in the question of determining the probability that a tuple of elements generate a 
finite group, which is the analogue problem where conjugacy is ignored, see for instance the 
paper |KLj of Kantor and Lubotzky. The paper [P] of C. Pomerance considers the question 
for abelian groups, when the conjugacy issue is also irrelevant, and his results do apply to 
our setting. The current paper will provide the beginning of the theoretical analysis of this 
type of algorithm for general finite groups. As a specific result, we will prove the following 
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result (informally stated; Theorem 16.11 gives the precise statement using the definitions of 
Section [2]): 

Theorem 1.1 (Boundedness of Chebotarev invariants for symmetric groups). There exists 
a constant c > with the following property: for all integers n ^ 1, the average number of 
independently, randomly chosen conjugacy classed of the symmetric group & n one must pick 
before ensuring that any tuple of elements taken from each of these classes generate & n , is 
at most c. In fact, for any k ^ 1, there exists c^ ^ such that the average of the k-th power 
of this number is bounded by Ck for all n. 

Here is the rough outline of this work: we consider probabilistic models in Section [21 and 
define an invariant, which we call the Chebotarev invariant of a finite group, using such a 
model (the name, based on the Chebotarev density theorem, is justified in Section [8]); it 
makes precise the informal notion in the statement of Theorem 11.11 Computing this invari- 
ant is seen to be related to very interesting questions of group theory, independently of any 
arithmetic motivation. In Section El we indicate how to compute this invariant for abelian 
groups (based on Pomerance's work) and in Section H] we consider solvable groups of a certain 
"extremal" type. In Sections and [71 we consider theoretical and numerical examples for 
non-abelian, often non-solvable, groups - in particular alternating and symmetric groups, 
proving Theorem 11.11 Finally, Section [8] makes some informal remarks concerning the appli- 
cability of our results to arithmetic situations (the original motivation); as we will see, there 
are non-trivial difficulties involved, and we hope to come back to these questions later. 

Notation. As usual, \X\ denotes the cardinality of a set, and F q is a field with q elements. 
If G is a finite group, and H C G, we write 

u G (H) = u(H) = ||f . 

We write G^ for the set of conjugacy classes of G, and for C C G", we also write u G (C) or 
v(C) for v(C), where C C G is the union of all conjugacy classes in C. 

In fact, as a matter of convenience, we will usually denote in the same way a subset of 
conjugacy classes and the corresponding set of elements in G, unless it is not clear in context 
if c G C means that c is a conjugacy class or an element of G (we will write often c" for a 
conjugacy class, avoiding most ambiguity). 

We recall that a geometric random variable X with parameter p G [0, 1] on a probability 
space is a random variable taking almost surely values in the set of positive integers, with 

(1.1) P{X = k) = p(l - pf- 1 
for k ^ 1. We then have 

(1.2) E(X)=p'£k(l- P ) k - 1 = -, E(X 2 ) = ^, VpO 



p p* p 2 



By / <C g for x G X, or / = 0(g) for x G X, where X is an arbitrary set on which / is 
defined, we mean synonymously that there exists a constant C ^ such that |/(x)| ^ Cg(x) 
for all x G X. The "implied constant" refers to any value of C for which this holds. It 
may depend on the set X, which is usually specified explicitly, or clearly determined by the 



This means distributed in proportion with the size of the conjugacy class. 
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context. Similarly, / xg means that / <C g and g <C / on the same set. On the other hand 
f(x) ~ g(x) as x — > xq means that f(x)/g(x) — > 1 as x — > xq. 

2. The Chebotarev invariant of a finite group 

In this section, we describe a natural probabilistic model for the recognition algorithm 
described previously. Fix a finite group G. We first remark that, whereas it does not make 
sense to say that a conjugacy class lies in a certain subgroup, unless the latter is a normal 
subgroup, it does make sense to say that it lies in a conjugacy class of subgroups. With that 
in mind, we define: 

Definition 2.1. Let G be a finite group, and let C = {Ci, . . . ,C m } C G' be a subset of 
conjugacy classes in G. Then C generates G if, for any choice of representatives gi G Ci for 
1 ^ i ^ m, the elements of the tuple (g±, . . . , g m ) generate G. Equivalently, C generates G if 
and only if there is no (proper) maximal subgroup H of G that has non-empty intersection 
with each of the Cj. @ 

The equivalence of the two definitions is quite clear contrapositively: if there are gi G Ci 
which generate a proper subgroup Hi, then each Cj intersects any maximal proper subgroup 
H of G that contains Hi, and conversely. Note also that the second condition can be stated 
by saying that there is a conjugacy class of maximal subgroups containing C. 

The following well-known lemma (due to Jordan) is the basic fact underlying the whole 
technique: 

Lemma 2.2. Let G be a finite group. Then the set G" of conjugacy classes generates H . 
In other words, there is no proper subgroup of G which contains a representative from each 
conjugacy class. 

Simple as this is, one should also recall at this point that the analogue of this lemma is 
false for infinite groups (even compact groups, such as SU(n), n ^ 2); for further discussion 
of various interpretations of this lemma, see [Slj . 

Now, let (f2, E, P) be a fixed probability space with a sequence X = (X n ) n ^i of G-valued 
random variables 

X n : i7 — y G, 

and let X\ be the conjugacy class of X n in G^: those are G"- valued random variables. 

Intuitively, those (X^) are the conjugacy classes that we see coming "one by one"; the 
Chebotarev invariant (s) looks at when we get enough information to conclude that those 
conjugacy classes can not all belong to some proper subgroup of G. 

We now define a random variable tx,g ( a waiting time) by 

(2.1) t x ,g = min{n ^ 1 | {X[, . . . generate G} G [1, +oo]. 

This depends on the sequence X = (X n ), and it may be always infinite (e.g., if X n = 1 
for all n\). But it is, in an intuitive sense, the "finest" invariant in terms of this probabilistic 
model. To obtain more compact and purely numerical invariants, it is natural to first take 
the expectation; this takes values in [1, +oo]. 

2 Alternately, following |Dlj . one says that elements (c/i, . . . ,g m ) invariably generate G if their conjugacy 
classes generate G in the sense above. 
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Definition 2.3. Let G be a finite group, X = (X n ) a sequence of G-valued random variables 
and tx,g the waiting time above. The Chebotarev invariant of G with respect to X, denoted 
c(G;X), is the expectation c(G;X) = E(tx,g) of this random variable. 

Remark 2.4. We focused on conjugacy classes because this is how applications to Galois 
theory are likely to arise, but of course one can similarly define an invariant using the 
original random elements (X n ) in G. If G is abelian, the two coincide. 

To have an unambiguously defined invariant, we must use a specific choice of sequence 
(X n ). The natural model is that of independent, uniformly distributed elements in G: if 
(X n ) are independent and identically uniformly distributed G-valued random variables, so 

and all n ^ 1, 
and all n ^ 1, 

then we call c(G;X) the Chebotarev invariant, and we just write c(G). 

Remark 2.5. It may be of interest, at least for numerical purposes, to use a sequence (X n ) 
which is not independent, but is obtained, for instance, by a rapidly mixing random walk 
on G. Also, the arithmetic analogue for computing Galois groups may be interpreted as 
involving non-independent and non-uniform choices of conjugacy classes (see Section [8]). 

Other numerical invariants may of course be derived from Tx,g, starting from the higher 
moments E(t xg ) for k ^ 1. Since it is probabilistically most important, when the expecta- 
tion of a random variable is known, to also have a control of its second moment, we define 
formally the secondary Chebotarev invariant: 

Definition 2.6. Let G be a finite group, X = (X n ) a sequence of G-valued random vari- 
ables, and let Tx,g be the waiting time above. The secondary Chebotarev invariant of G with 
respect to X, is the second moment c 2 (G;X) = E(r|- G ). If (X n ) is a sequence of indepen- 
dent, uniformly distributed random variables, then we write 02(G) and call it the secondary 
Chebotarev invariant. 

We will now give formulas for the two Chebotarev invariants (in the independent case), 
which are expressed purely in terms of group-theoretic information. This is useful for explicit 
computations, at least for groups which are very well understood (but often the probabilistic 
origin of c(G) should also be kept in mind.) 

To state the formulas, we must introduce the following data and notation about G. Let 
max(G) be the set of conjugacy classes of (proper) maximal subgroups of G (if G is trivial, 
this is empty); for a conjugacy class of of maximal subgroups IK G max(G), let K" denote 
the set of conjugacy classes C of G which "occur in IK", i.e., such that C fl Hi 7^ for some 
Hi in the conjugacy class kJI Moreover, if I C max(G) is a set of conjugacy classes of 
maximal subgroups, we let 

(2.2) ^5=0 ^ 



uxia-i 

P(X n = g) = jir for all g G G, 

and hence ^ 

P(Xl = g*) = j|l for all g* G G\ 



Note that this depends on the underlying group G. 
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the set of conjugacy classes of G which appear in all subgroups in /. 
Then we have: 



Proposition 2.7. Let G be a non-trivial finite group. With notation as above, we have 

i 

[)' 



(2.3) c{G)= Y 1 ! „ 



JCmax(G) 



and 



/Cmax(G) V 11 K IJ /Cmax(G) V V 1,1 

These formulas do not apply for the trivial group, since they lead to empty sums which 
are zero, whereas the definition leads tcQ c(l) = 1, c-i{X) = 1. 

Probabilists will have noticed that the first formula (at least) is very similar to the one for 
the expectation of the waiting time for a general coupon collector problem. There is indeed 
a link, which is provided by the next lemma where independence of the random elements X n 
is not required. 

Lemma 2.8. Let G be a non-trivial finite group and X = (X n ) a sequence of G-valued 
random variables. The waiting time tx,g ^ s equal to 

t x ,g = max f K , 

IK £ max G 

where 

(2.5) foe = mm{n ^\\X\i M»} ; 

note that f<y{ depends also on X. 

In other words, tx,g is a l so the maximal n such that we need to look at Xi for i up to 
n, before we witness, for every conjugacy class "K of maximal subgroups, some X n which is 
incompatible with the groups in this class 'K. This is very close to a coupon collector problem 
(see, e.g., |FGT] for a general description of this type of problems), where the "coupons" we 
need to collect correspond to the conjugacy classes which are not in CK", as 'K ranges over 
max(G). But since a single X n may serve as coupon for more than one IK", this does not 
exactly correspond to standard coupon collector problems!! Because of this, we state and 
prove the following general abstract result, which may have other applications. 

Proposition 2.9. Let (Q, E, P) be a probability space, D a finite set. Let (Z n ) be a sequence 
of D -valued random variables. Let £ be a non-empty finite collection of non-empty subsets 
of D, and let 

f£ = min{n 1 | for all E G £, there exists some m ^ n with Z m G E} 



4 One may argue that there is no need to look at any elements to be sure of generating the trivial group, 
but this does not correspond to the definition. 

5 This has been called the "coupon subset collection problem" by Adler and Ross [X] , 
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be the waiting time before all subsets E G £ have been witnessed in the sequence (Z n ). For 
I C £ ; non-empty, let 

(2.6) Tj = min{n ^ 1 | Z n E E for some subset E e J}. 

(1) Assume Tj < +oo almost surely for all non-empty subsets I C £. TTien we nave 

(2.7) r £ = 

0^/c£ 

(2) Assume the Z n are independent and identically distributed random variables and let \i 
be their common law. We have 

f2Sl Ef^-V ( " 1)|J ' +1 -V 

and iae subsets in £ are disjoint, we have 

»+oo 



(l-IJa-expC-M^)*)))^. 

£e£ 

(3) Pfe have 



£ /i(U Se /^)V A*(U B6 /J5) 

7^0 

When £ is the set of singletons in .D, where we have exactly the coupon collector problem, 
the formulas for the expectation are well-known (see, e.g., |FGTl Theorem 4.1] for the integral 
formula); we have not seen general formulas for the second moment in the literature. 

Proof of Proposition \2.9\ . To simplify notation, define 

(2.11) E T = |J E, 

Eel 

for each I C £. Formula (12. 7p - which implies in particular that T£ is finite almost surely - 
can be checked, e.g., by seeing T£ as the length of the subset 

\J[0,T {E} ]cR, 

£e£ 

(which is therefore almost surely finite by assumption on the Tj), and applying the inclusion- 
exclusion formula for the measure of a union of finitely many sets (for the Lebesgue measure, 
or the counting measure on Z, equivalent ly): 

|U[°' T w]|= E c— p| [o,^^>]|, 

£e£ 0^/c£ Eel 

at which point it is enough to observe that, for any / C £, we have 

| p| [0,T {e} ]| =mmT {E} = T 7 . 
Eel 

We can now finish the computation of E(rg) in (2), in the case of independent random 
variables. Indeed, in that case the random variable Tj is distributed like a geometric random 
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variable with parameter p = P(Z n G Ej) (see (11. ip ) for any non-empty subset / C £, so 
that taking expectation in (12. 7p and applying (II. 2p . we obtain the result. 

The integral expression (I2.9p is a consequence of (12. 8p and the additivity of measure for 
disjoint sets: it suffices to expand the product and use 



-oo 



e at dt = -, for a > 0. 
a 

Finally, to compute the second moment in the independent case, we start with the same 
formula (12.71) to get 

E (^ 2 ) = EE(- 1 ) |/|+|JlE ( T ^)- 

0^/C£ 
0^Jc£ 

We first transform this by applying (I2.16P in Lemma r2.14l below to compute E(T/Tj). This 
gives 



* = \^Jzi»_>_ + _L_ 



E(r|) = EE 



0^/c£ 
0^Jc£ 



( 2 ' 12 > (by symmetry). 



07^/C£ 
0^Jc£ 

To continue, consider more generally arbitrary complex coefficients /?(/) defined for / C £, 
and the expression 

f_i )W+\J\ 

0^Jc£ 

note that E(r|) is a simple combination of two such expressions. 

We proceed to reduce W(f3) to a single sum over I C £ by rearranging the sum according 
to the value of / U J: 

£ ^EEn" «')■ 

IU.J=K 



The inner sum is rearranged in turn as 

^^(-l)IWI/3(/)= £ (-l) |J| /3(/) E ("^ 

0^/,JC£ 0^/Cif 0^JCK 

IUJ=K IUJ=K 



= E (-i) m+l Mw E (- 1 ) 17 ' 1 

r'u(if-/)^0 

since the subsets J with J U J — K are parametrized by I' C / using the correspondence 
J' H> (A' - J) U J' with inverse J ^ J Hi. 

For fixed J, the last summation condition /' U (K — I) ^ is always valid, unless I = K. 
In that last case, it only excludes the set I' = from all I' C /. Since we have, for any finite 
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set X, the binomial relation 

£(-i)' y| = o, 

YcX 

it follows that the double sum is simply given by 

x;e(- i ) |/i+iji ^ / ) = (- i ) |jfi+ M. 

IUJ=K 

and hence 

Applied to the expression (12.121) . this leads precisely to ( 12 . 1 Of) . □ 



To deduce Proposition 12.71 we apply this proposition with 

Z n = Xl D = G i £ = - JC tt | IK G max(G)}, 

in the case where the (X n ) are independent uniformly distributed on G, so that the common 
distribution is fi = v. Since, for I C maxG, we have 

■Kei Me/ 
the formulas (12. 8p and (12.101) give exactly the claimed formulas (12.31) and (12 .4p . 

Remark 2.10. Note the following strange-looking "linearity" property, which can be checked 
from our formulas and (12. 16ft : for 1, we have 

c(G) = E(r G )= Yl (-l)' /|+1 Em), 

0^/CmaxG 

c 2 (G) = E(r G ) = Yl (-1) |/|+1 E(T 7 2 ). 

Remark 2.11. These formulas can only be useful for practical computation if the number 
of conjugacy classes of maximal subgroups of G is fairly small, or if they are very well 
understood. As a theoretical instrument, they suffer from the fact that it is very hard to 
use them to guess or estimate the actual value of c(G). For instance, it is not clear how to 
recover even the trivial lower bound 

(2.13) c(G) > 5(G), 

where 6(G) is the minimal cardinality of a generating set of G. We will give examples later 
on where this bound is very close to the truth, even with groups of size growing to infinity. 

Remark 2.12. Another natural formula is 

(2.14) c(G) = l + J2 P ( X i 8 ' • • • ' X " do not £ enerate G )> 

which is also valid for G — 1. 

Indeed, since tq takes positive integer values, we have the familiar formula 

E(r G ) = YP(r G >n), 
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and clearly 

{t g ^ n} = {X\, . . . , Xl_ x do not generate G}, 

for n ^ 1. When n = 1, this is the certain event, with probability one, thus leading to (I2.14p . 

A moment's thought shows that one can also identify this formula with the one coming 
from (12. 3p by expanding the geometric series: 

E i _ JT'w) = E (-i)'" +1 E"(n*"r 

/Cmax(G) Vl l:Ke/ ^ /Cmax(G) ™>0 JTG/ 

Ii4 



= 1 + £( £ (-i) |/|+ vn^)' 

™>1 /Cmax(G) 

(where the term with n = is only equal to 1 for max(G) ^ 0, i.e., G ^ 1). 
This is not really a different proof of Proposition 12.71 since the relation 

P(Xl,...,Xl do not generate G) = f| 

/Cmax(G) Hel 



is proved by inclusion-exclusion, exactly as in the proof of Proposition | 

The point of (12. 14ft is rather that it leads to another lower bound (for 1): 



c(G) + p ( x l • • • , X* do not generate G) 

for any fixed k ^ 1, and this may be quite useful because there has been a large amount of 
work on the estimation of those probabilities when k is small, e.g., k = 2 if G is not cyclic. 
For instance, Dixon (for G = A n and n — > +oo) and Kantor and Lubotzky (for G = G(F 9 ) 
a simple classical group and q — > +oo) have shown that in those cases we have 

P(Xi,X% do not generate G) — > 

as n (resp. q) goes to infinity, indeed with quantitative estimates (see |KL] ) - but note the 
probabilities with elements and with conjugacy classes may behave rather differently (e.g., 
for G = PSL(2, F p ), the probability that two random elements generate G is very close to 1 
for large p, but there is a probability converging to 1/2 that two random conjugacy classes do 
not generate G, see the proof of Theorem 15.11) . In a sense, the Chebotarev invariant is thus 
a refinement of these type of probabilities. We refer to |D2] for a brief survey of probabilistic 
Galois theory, and to |D1] for the analysis of the case of symmetric groups. 

Remark 2.13. As explained in |S1[ Th. 5], we have 
(2.15) 

for any conjugacy class of maximal subgroup of G (this is due to Cameron and Cohen). 

We now prove the lemma which supplies the formula (12. 16ft used in the proof of the 
proposition. 



Lemma 2.14. With notation as in Proposition ^.^ in particular with independent identically 
distributed random variables (Z n ), for any two non-empty subsets I , J in £ ; we have 

( 2 - 16 ) HTiTj) = — J—f—L- + 



fi(E IUJ ) \fi(Ei) fi(Ej) 

Proof. This is a fairly direct computation, but not very enlightening (at least in our presen- 
tation; there might be other approaches that makes this more transparent). 

We first compute the joint distribution of Tj and Tj, and for this, we use the shorthand 
notation 

P = K E i)> Q = K e j)> r = n(E InJ ), s = fi(E IUJ ), 
p' = fi(Ej-Ej), q' = /i(Ej — Ei). 

Then we have (generalizing the geometric distribution of a single Tj): 

{(1 — s) n ~ 1 r, if n = m ^ 1, 

(1 - s) m -Y(l - p)"^ x p, if n > m ^ 1, 
(1 - s) n -y(l - q^-^q, ifm>ra^l. 

To justify, e.g., the second of these, note that Tj = n > m = Tj means that must 
not be in E^jj for k ^ m — 1, Z m must be in Ej but not in Ej, Z^ must not be in Ej 
for m < k < n, and finally Z n must be in Ej\ then the independence of the (Z n ) gives the 
formula. 

Now we write 

E(T/Tj) = nmP(Ti = n and Tj = m), 

and we split the sum according to the three cases, say 

E(T I Tj)=Q 1 + Q 2 + Q 3 . 
Introducing further the functions 



x) n -\ 



we have the expressions 



Qi = J2 n2p ( T i =Tj = n)=rY^ n 2 (l 



8)^ = ^2(8) 



and 



Q 2 = nmPiTj = n and Tj = m) 

l^m<n 

= q 'p Mi-sr'ti-pr" 1 " 1 

= q'p H 1 - s ) m_1 ^ (m + A;)(l - p) k ~ l 



= q'p(Mp) 1 P2(s) + Vi(p)yi(s)), 

while Qs is given by the same expression after exchanging p and q, p' and q' . 
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Since, by Taylor expansion, we have 



we obtain 



1 1 2-x 2 1 

<P0{X) = - tfl{x) = "J, = — — = -Z~-2- 

JU JU tXj JU ikj 



E( r r T J ) = i(? + ^ + ^-i, : 

s Vs ps qs 



by adding the three terms. Finally, the relations 

q' + p = s, p' + q — s, 

lead to the simplified expression 

^/r^^x 1/2 q' p' \ 1/2 s—p s — q 
E(TjTj) = _(_ + ^- + £--l) = _(_ + 1 + 1 _ l 

s Vs ps qs / s Vs gs 

1/21111 
= -(- + + 1 

s Vs p s q s 

which gives ( I2.16p . □ 

We now present some easy formal properties of the Chebotarev invariants (attached, unless 
stated otherwise, with a sequence of independent uniformly distributed random variables). 

The first lemma may be used to simplify and expand the range of groups covered by certain 
computations (this is also observed by Pomerance [P] , for the case of numbers of generators 
instead of conjugacy classes): 

Lemma 2.15. Let G be a finite group, and let $(G) be the Frattini subgroup ofG, i.e., the 
intersection of all maximal subgroups of G. Then, for any normal subgroup N <G such that 
N C in particular for N = $(G), we have 

c(G) = c(G/N), c 2 (G) = c 2 (G/N). 

Proof. Let H = G/N. Moreover, we have $(H) = $(G)/N and hence H/<$>(H) ~ 

This means that we need only prove the result when N = $(G), the general case following 

by applying this to H. 

Let 7r : G — > G/$(G) be the quotient map. If (X n ) is a sequence of independent random 
variables uniformly distributed on G, the Y n = n(X n ) are independent and uniformly dis- 
tributed on G/<f>(G). Moreover, for any n ^ 1, the elements (X[, . . . ,X%) generate G if and 
only if the elements (Yf, . . . , Y^) generate G/&(G): indeed, this follows from the basic fact 
that a subset S C G generates G if and only if tt(S) generates G/$(G) (this is applied to 
all sets S = {xi, . . . , x n } where X{ conjugate to Xi). This gives the result immediately from 
the definition of the waiting times. □ 

We next consider products: 

Proposition 2.16. Let G±, G 2 be finite groups such that the only subgroup H C G\ x G 2 
which surjects by projection to both factors is H = G. Then we have 



c{Gr xG 2 )^c{G 1 ) + c{G 2 )-l. 
n 



Examples of groups G\, G 2 satisfying the hypothesis are any pair of non-isomorphic simple 
groups; note that this proposition suggests that sometimes c'(G) = c{G) — 1 would be a more 
natural invariant to consider, since we then have the simpler inequality 

c'(G 1 x G 2 ) ^c'(G 1 ) + c'(G 2 ). 

Proof. With G = G\ x G 2 and denoting X n = (Y n , Z n ) 6 Gi x G 2 a sequence of independent 
uniformly distributed random variables, it is clear that (Y n ), (Z n ) are similarly independent 
uniformly distributed on G\ and G 2 respectively. We then have the inequality 

t g ^ max(r!, r 2 ) ^ n + r 2 - 1 

(since Tj ^ 1 and max(m, n) ^ n + m — 1 for integers n, m ^ 1), with 

T\ = min{?7, ^ 1 : (Yf, . . . , Y^) generate Gi}, r 2 = mm{n ^ 1 : (Z[, . . . , Z^) generate G 2 } 

which are distributed like tq x , t"g 2 (indeed, if n ^ max(ri,r 2 ), then the group generated 
by any elements in X| = (Y£, Z|) surjects to G\ and G 2 , hence it must be equal to G by 
assumption). Taking expectation, we get the inequality stated. □ 

The next result gives upper and lower estimates for the Chebotarev invariant using smaller 
sets of maximal subgroups than max(G). This can be very useful in particular for the 
asymptotic study of c{G n ) for a sequence of finite groups (G n ), as we will see later on. 

Proposition 2.17. Let G be a finite group, and let M C max(G) be an arbitrary non-empty 
finite subset of maximal subgroups. Let 

f M = maxf M . 

JteAf 

with notation as in (12.51) and 

(2.18) p M = „(Gt- |J K» 

JCemax(G)-M 



We then have 



and 



E(f^) ^ c 2 (G) < E(t* ) + - 1. 

Pm 

Proof. Define the additional waiting time 

r* = min{n ^ 1 | X n £ (J 

We then note the inequalities 

?m < r G ^ max(f M , r*) ^ f M + r* - 1, 

where the first inequality is obvious, while the second follows because, for n = max(fM> t*), 
we know that the group generated by (X{, . . . ,X%) is not contained in any subgroup in a 
conjugacy class of maximal subgroups CK e M, and that this group also contains one element 
which is not conjugate to any element in a subgroup not in M. 
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Now we take expectation on both sides; observing that, by independence, r* is distributed 
like a geometric random variable with parameter p m given by (I2.18p . we obtain the first 
inequalities, using Proposition 12.91 and (jl.2p . 

Similarly, for the secondary invariant, we use the inequalities 



and get 



f 2 M «C r 2 <k max(f M , t*) 2 ^ f 2 M + (r*) 2 - 1, 



E(f A 2 f ) < c 2 (GQ < E(f2r) + E((r*) 2 ) - 1 = E(f 2 / ) + 1™ - 1. 

Pa/ 



□ 



The following immediately follows: 



Corollary 2.18. Let (G n ) be a sequence of non-trivial finite groups, and let v n denote the 
corresponding density. For each n ^ 1, let M n be a non-empty subset ofm&x(G n ), and 
assume that 

(2.19) lim V J II '}<f\ = 0, 

MGmax(G„)-A/„ 

i.e., the proportion of elements represented by a conjugacy class in some subgroup in M n 
goes to zero. Then we have 

c(G n ) = E(f M J + o(l), and c 2 (G n ) = E(f 2 M J + o(l), 



as n — >■ +oo, with notation as in Proposition 2.11 



The following sections will now take up the problem of computing, or estimating, the 
Chebotarev invariants for various classes of groups. 

3. Abelian and nilpotent groups 

In this section, we look at finite abelian and nilpotent groups G. In fact, because nilpotent 
groups have the (characteristic) property that [G, G] C $(G) (see, e.g., \Ro\ Th. 11.3, (v)]), 
Lemma [2. 151 shows that if G is a nilpotent group, we have 

c(G) = c(G/[G, G\), c 2 (G) = c 2 (G/[G, G]) 

which are Chebotarev and secondary Chebotarev invariants of abelian groups. This applies, 
in particular, to all p-groups. 

We will not use the formula from Proposition 12.71 (although it is possible, as was done in 
a first draft, to do some computations using it), because in abelian groups there tends to 
be many maximal subgroups up to conjugacy - since conjugacy is now trivial. Following 
the work of Pomerance [P], who computed c(G) (with different terminology) for any abelian 
group G, we will use another description of the Chebotarev waiting time in the case of 
abelian groups. 

Theorem 3.1 (Pomerance). Let G be a finite abelian group, and for any prime number 
p | let r p (G) = dim.F p (G/pG) be the p-rank of G. Let 5(G) = maxr p (G) be the minimal 
cardinality of a generating set ofG. Then we have 

c(G)=5(G) + J2(l-U II (l-^ (G)+ - l) ))- 

j>1 p||G|l^r- p (C?) 
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In particular, for G = Z/nZ with n ^ 2, we have 

Hid) 



(3.i) c ( G ) = -Ei 



al- 1 

d\n 

and for G = F*, where F p = Z/pZ, with p prime and k ^ 1, we have 
(3.2) -^T- 

This is [Pj Theorem] and immediate corollaries of it. The formula for G = Z/nZ might be 
easier to get directly from Proposition 12.71 Indeed, the subgroups of Z/nZ are the groups 
H d = dZ/nZ for d \ n, with v(H d ) = d~ l and H d r\ H e = H^e], an d the maximal subgroups 
among these correspond to minimal divisors of n for divisibility, i.e., to the primes p dividing 
n. Then a non-empty subset / of max(G) can be parametrized by the corresponding subset 
of prime divisors of n, or equivalently by the squarefree divisor d > 1 of n which is the 
product of those primes. In this correspondence, we have 

f)H = f)H p = H d , hence u ( f| if) = u(H d ) = 



hgi p \d Hei 



and (-1)I J I = fj,(d), hence O gives the stated formula for c(Z/nZ). 

We have similar results for the secondary Chebotarev invariant; Pomerance mentions the 
possibility of computing these, but does not give any results in his paper. 

Theorem 3.2. Let G be a finite abelian group. With notation as in Theorem \3.1[ we have 
c 2 (G) = 5{Gf + £ (2j + 25(G) - 1) (l - J] J] (1 - p-W^-O)) . 

J>1 p\\G\ l<i<r p (G) 

In particular, we have 

C2 (Z/nZ) = -5>(rf)7T^^ 

for n ^ 1 anci 

c 2 (Fj) = c(Fj) 2 + £ 



(pj - 1)2' 



/or p prime and k ^ 1 . 



Proof. The first result is obtained by reasoning as in (Pj p. 195], with r and (r + j) there 
replaced by r 2 and (r + j) 2 . The point is that he shows that 



P((X h X s{G)+j ) generate G) = ] (1 - p 



{8(G)-r p {G)+j+iY 



p\\G\ lC^rp(G) 
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To deduce the values for G = Z/p k Z, it is simpler to use the description^ 

k 

where the Gj are independent geometric random variables with parameters pj = 1 — p~K 
Concretely, they can be defined as follows 

G k = min{n ^ 1 | X n ^ 0}, 
= min{72 ^ 1 | dim Fp (X Gk+n , X Gk ) = 2}, ... 
Gi = min{ra ^ 1 | dim Fp (X G2+n , X G2 , . . . , X Gfe ) = fc}, 

which, by independence of the (X n ), are easily checked to be indeed independent geometric 
variables with the stated parameters. 

This decomposition leads to the formula for 02(G) immediately, using (11.21) and additivity 
of the variance of independent random variables. □ 

The formula of Pomerance gives a quick way to understand the limit values of Chebotarev 
invariants for abelian groups with a given rank 5(G). 

Corollary 3.3 (Pomerance). For any fixed integer k ^ 1, and any abelian finite group G 
with 5(G) = k, we have 




In particular, the Chebotarev invariants for cyclic groups are bounded. 
Corollary 3.4. For any fixed k, we have 

c(F k p ) = k + 0(p~ 1 ), c 2 (Fj) = k 2 + 0(p~ 1 ), 

and 

P(r n ^k)<€.p~\ 
where the implied constants depend only on k. 

This last result shows that, for vector spaces over a finite field, the Chebotarev invariant 
is strongly peaked around the average, which is itself close to the dimension. 

Proof. Only the last inequality needs (maybe) a bit of explanation. Since r F k takes positive 
integer values ^ k, we have 

|r F fc -k\^l 

if r F fc 7^ k. Hence, if r F k 7^ k, we have 

\r n - c(F k p )\ > \r F k -k\- \c(F k p ) - k\ > 1 - \c(F k ) - k\, 
and if we furthermore we have p ^ po, where po (depending on k) is chosen so that 

k < c(F k ) <:k + 1/2 

6 Which is the analogue of the decomposition of the waiting time for the standard Coupon Collecting 
Problem in a sum of geometric random variables. 
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for all p ^ po, it follows that 

{r n ^ k} G {\r n - c(F k p )\ ^ 1/2} 

for such p, and then the Chebychev inequality gives 

P(r F ,^A;)^4V(r F ,)«^ 1 

for p ^ po, where the implied constant depends on k. Increasing this constant if needed 
(e.g., taking it to be at least po), we can also claim that this inequality holds for p ^ 2. □ 

Remark 3.5. In particular, for cyclic groups, the Chebotarev invariant is at most, and its 
limsup is, the constant 



(l - tttt) = 2.705211140105367764. 



This asymptotic behavior is not without interest (and some surprise): on the one hand, 
we see that c(Z/nZ) remains absolutely bounded, despite the existence of cyclic groups with 
many subgroups, and on the other hand, we see that it is not always close to the minimal 
number of generators. 

Concerning the first point, notice for instance that a "naive" invariant is given by 

E •w-'Ll 

ff£max(G) p\n 

and this is unbounded as n grows for G = Z/nZ (though it is <C log log log n, and thus 
bounded in practice...), see the discussion surrounding (19 .2p below for occurrences of such 
quantities instead of the Chebotarev invariant. 

For the second, note that (interpreting 1/C(1) = 0) the limsup we found is also N where 

^ = D (i - C(A0 -1 ) 

is sometimes called the Niven constant. Niven obtained it as the mean-value of the maximal 
exponent of a prime dividing a positive integer: 



:rr E °u> 



N = lim 

n— >+oo n 

with 

a(j) = max{u ^ | p v \ j for some prime p} 

for j ^ 2 (see [N]). The explanation for this coincidence is that ((k)" 1 , for k ^ 2, is both 
the (asymptotic) density of primitive vectors in Z k and that of fc-power-free integers. 

Remark 3.6. If one uses Proposition 12.71 instead, one can prove (after some computation) 
the following formulas 

P 



S l-p-3 \j) P ' 

^)=E(-^(^CV' (I " 1)/2 
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where 

'k\ (1 -/)••• (1 - p k ~ j+1 ) 



JJ P (l-^')---(l-p) 
are the p-binomial coefficients. Note that those formulas do not immediately reveal the lim- 
iting behavior as p — > +00, since the summands have different degrees as rational functions 
of p. 

4. A SOLVABLE EXAMPLE 

The results of the previous section, as well as those we will see in the next one, reveal (or 
suggest) rather small values of the Chebotarev invariants, in comparison with the size of the 
groups. The following example in the solvable case exhibits very different behavior (possibly 
the worse possible). 

Proposition 4.1. For q a power of a prime, let 

it U 



q \\0 1 I 1 ? ' q 

be the group of translations and dilations of the affine plane ~F 2 q of order q(q— 1), isomorphic 
to a semi-direct product F 9 >iF ? x . 
(1) We have 



(4.1) c(H g )=q-q- 1 £ - 



fi(d) 



(4.2) c 2 (H q ) = q(2q-l) + c 2 (Z/(q-l)Z) + £ ^(d) * + 

(2) For q ^ 2, we have 

(4.3) c(if ff ) = q + 0(r(g - 1)), c 2 {H q ) = q(2q - 1) + 0(r(g - 1)), 

where r(q — 1) is £/ie number of positive divisors of q — 1. In particular, c(H q ) ~ g os 
g — >■ +00. 

Since we have a split exact sequence 

1 _> f, -> F * _> J 

and the two surrounding groups are isomorphic to F*, where q = p k with p prime, and 
to a cyclic group Z/(g — 1)Z, with Chebotarev invariants tending to k as p gets large, and 
bounded, respectively, this shows in particular that the Chebotarev invariant can jump quite 
uncontrollably under extensions. 

The proof will use Proposition 12.71 We start with a lemma that is certainly well-known, 
but for which we give a proof for completeness and lack of a suitable reference. 

Lemma 4.2. (1) There are q conjugacy classes in H q ; they are given, with representatives 
of them, by 

9b =(oi)' 9b = {9^H q \ det(g) = 6}, \g\\ = q, 
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where & G F 9 X - {1}, and 

ld=(l fj, Id« = {Id}, 

«=(J}J, = {g eH q - {Id} | det(g) = 1}, | U »| = g-l. 
(2) The conjugcay classes of maximal subgroups of H q have representatives given by 



I 1 J I aGF « 



= U a n UeH a \ ae(F*Y andte 



and 

where i runs over the prime divisors of q — 1 . 
Proof. (1) We have the general conjugation formula 



(4.4) 



a t\ ( b v\ fa t\ 1 _ I b av + t(l — b) 



1 10 1 10 1 10 1 



from which it immediately follows that all elements with det(g) = 6^1 are conjugate, and 
gives therefore the q — 2 conjugacy classes with representatives gb described. For b = 1, it is 
clear that all elements with b — 1, v ^ form the conjugacy class it", and only the identity 
class remains. 
(2) Denote 



U = H q nSL{2,F q ) = {(l fj | t6F,}. 



Let H <Z H q be a maximal subgroup. Let D = det(H) be the image of the determinant 
restricted to H. Since 

H C det^p), 

we have either det" 1 (-D) = H q , i.e., D = F*, or if = det _1 (JJ). In this second case, the 
subgroup D must be a maximal subgroup of F^ for H to be maximal, which implies that H 
is of the form Ci for some £. Conversely, such a subgroup is maximal because if we add any 
extra element g and let H' = (Ci,g), the fact that U C C% implies that some 

a' 
1 

with a' (jz (F* ) e , is in H', and then by maximality in F£ , we have H' = H q . 

Note that the Ce are normal in H q , hence also pairwise non-conjugate. 

In the first case, when D = F^ , i.e., when the determinant restricted to H is surjective, we 
claim that the determinant is also injective on H: indeed, otherwise, there exists u G U DH, 
u ^ 1, say 

u — \ q * J with t 0. 
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For any a G F* by surjectivity there exists a (a) G F q with 



a a(a) 
1 



and by applying the relation 



a a(a)\ (\ t\ fa a(a)\ _ fl at 
1 ) ' [o 1 J ' \0 1 J = [o 1 

for all a G F* we conclude that in fact U C H. Then \H\ is divisible both by q and by 
q — 1, hence H = H q , contradicting the assumption that if is a proper subgroup of H q . So, 
in this second case, the determinant gives an isomorphism H ~ F^ . Then a generator of H 
is in one of the conjugacy classes (it has distinct eigenvalues in F£ ), hence it is conjugate 
to an element (generator) of A, and H itself is therefore conjugate to A. □ 

Proof of Proposition \4-l\ First of all, in addition to the maximal subgroups Ci given by 
Lemma 14.21 there are subgroups Cd for all squarefree divisors d \ q — 1, the inverse image 
under the determinant of the subgroup of order (q — l)/d in the cyclic group F*. 
Given a subset I C ma.x(H q ), we now compute the density of conjugacy classes in 



w} = n ^> 

as follows: 

- If A G I, then with V = 1 — {A}, and d the product of those primes I for which Cg G I' 
(including d — 1 when i 7 = 0), we have 

z/(!Kj) = — — g^ 1 and in particular = 1 — g" 1 . 

Indeed, we have to find the density of those elements of H q which are diagonalizable with 
eigenvalues 1 and a G Cd- These are exactly the conjugacy classes g\ with b G Cd — {1}, and 
the trivial class, so 

^ = l + ((g-l)/rf-l)g = g(g-l)/rf-(g-l) = 1 _ 1 

7 g(g-i) d q 

- If A ^ 7, then / corresponds to a divisor d \ q — 1, d 1, and we have 

"(ft}) = ^ 

since we must now compute the density of elements of H q which have det(g) G Cd, and this 
is 

g((g-l)/rf-l) + l + _ 1 

g(g-l) d' 

Applying (12. 3p . and isolating the contribution of I = {A}, leads exactly to (14.11) and 
to dOJ). 

To deduce (14. 3[) for c(H q ), we may assume q = p k with p an odd prime, since for q even, 
we have 

c(H q ) = q + c(Z/(q-l)Z)=q + 0(l) 
19 



by Corollary 13.41 So for q odd, we write 



c(H q ) = q + c(Z/(q - 1)Z) - A(g) = q - A(q) + 0(1) 

where 

A(q)= y ^ . 

Since 1 — <i -1 + g -1 ^ 1 — > 0, we can bound this from above by 



and then proceeding as in the proof of Corollary 13.41 we obtain 

k^O p\q-l 

= 0(r(g-l)) 

since the series converges absolutely again. 

Finally, the asymptotics for C2(H q ) are obtained by essentially identical arguments. □ 

The proof confirms the intuitive fact that the large size of c(H q ) is due directly to the 
existence of a fairly small diagonal subgroup A (of index q) that contains elements conjugate 
to a very large proportion of elements of H q . So the waiting time is quite close to the waiting 
time until a non-diagonalizable element is obtained, which is a geometric random variable T 
with 

1/ l\ k ~i 
P(T = k) = - ( 1 - - J , for k ^ 1 
q V g/ 

(since very often, it will be the case that sufficiently many diagonalizable elements will have 
appeared by the time an element of U appears to generate the whole group). 

This is confirmed by the large second moment C2(H q ): it corresponds to a standard devi- 
ation of the waiting time which is 



'c 2 (H q ) - c(H q ) 2 ~ q, as q -> +oo, 
i.e., very close to the expectation, similar to the fact that 



V(T) = ^l-i. 

The groups G = H q also show that the inequality (12 . 151) is best possible (with the maximal 
subgroup H = A), as observed also in |Slj . so it is not surprising that they lead to high 
Chebotarev invariants. 
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5. Some finite groups of Lie type 



For specific complicated non-abelian groups, the Chebotarev invariant may be hard to 
compute exactly, except numerically using the formulas of Proposition 12.71 when feasible 
(we will give examples from computer calculations in Section [7]). However, if we consider 
infinite families of non-abelian groups, it may be that the subgroup structure is sufficiently 
well-known, simple and regular, that one can derive asymptotic information. In fact, using 
results like Proposition 12.171 it is not needed for this purpose to have complete control over 
all maximal subgroups. 

We illustrate this first simplest family of simple groups of Lie type. 

Theorem 5.1. (1) We have 

c(PSL(2,F. p )) = 3 + 0(p- 1 ), c 2 (PSL(2,F p )) = ll + O^" 1 ), 

for primes p ^ 2. 

(2) For all k ^ 2, we have 

P(r PSL (2,F p ) = k) = + Oip- 1 ) 

where the implied constant depends on k. 

(3) The same results hold for SL(2, F p ), and in fact 

(5.1) c(SL(2,F p )) = c(PSL(2,F p )), c 2 (SL(2, F p )) = c 2 (PSL(2, F p )) 

for all p. 

Note that the limit of c(SL(2, F p )) is not the minimal number of generators of SL(2, F p ) 
(which is 2, since SL(2, F p ) is generated by the two elementary matrices with 1 over and 
under the main diagonal.) 

For the proof, we will not use the formula of Proposition 12.71 although this could be done 
at least to prove (1) (since the subgroups of PSL{2, F p ) are well understood since Dickson, 
see, e.g., |Gil Th. 2.2]). Instead, we use the following criterion of Serre |S31 Prop. 19] (which 
is itself based on knowing the subgroup structure) . 

Lemma 5.2 (Serre). Let p ^ 5 be a prime number. Assume that G C SX(2, F p ) is a 
subgroup such that 

(1) The group G contains an element s such that Tr(s) 2 — 4 is a non-zero square in F p , 
and such that Tr(s) ^ 0; 

(2) The group G contains an element s such that Tr(s) 2 — 4 is not a square in F p , and 
such that Tr(s) ^ 0; 

(3) The group G contains an element s such that Tr(s) 2 G F p is not in {0, 1,2,4}, and is 
not a root of X 2 - 3X + 1. 

Then we have G = SL(2, F p ). 

Proof of Theorem \5.1[ We first notice that we need only check (15.11) and then consider the 
case of SX(2,F P ). These equalities are consequences of Lemma 12.151 sinc^ {±/} Q $ P , 
where <3> p is the Frattini subgroup of SX(2,F p ). Indeed, we may of course assume that 



In fact, it is known that there is equality, but we do not need this stronger fact. 
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p 2; then, if p is such that —I ^ $ p , there exists a maximal subgroup H of SL(2, F p ) 
which surjects to PSX(2, F p ). We would then have 

5L(2,F p ) = {±7} x 77 

which is impossible, since SL(2, F p ) is generated by the elements 

oO (1? 

both of which are of odd order p, hence contained in H (compare with jS2| IV-23]). 

Now we consider SL(2, F p ), and we assume p ^ 5. Let r = tsl(2,f p ) denote the cor- 
responding waiting time, and let n, T2, T3 denote the waiting times for conjugacy classes 
satisfying the conditions (1), (2) and (3) in Lemma [5.2[ i.e., for instance 

n = min{n > 1 : s = X s n has Tr(s) ^ and Tr(s) 2 - 4 is in (F* ) 2 }. 

Let also rj", be the waiting times for conditions (1) and (2) without the condition 
Tr(s) 7^ 0. Note that (1) and (2) are exclusive conditions. Moreover, each 7$ is a geometric 
random variable, with parameters, respectively 

(5.2) Pi = l + 0(p~ 1 ), p^^ + Oip- 1 ), p 3 =l + 0(p- 1 ), 

and for r*, r|, the parameters are also 

A = l + 0(p- 1 ), p; = l + 0(p- 1 ); 

all these facts can be checked easily, e.g., by looking at tables of conjugacy classes in SL(2, F p ) 
(for instance, |FH] p. 71]). 
We then have 

max(r*,r*) ^ r p < max(ri, r 2 , r 3 ), 

where the right-hand inequality comes from Lemma 15.21 and the left-hand inequality is due 
to the fact that the Borel subgroup 



M(o r')} CJ W 



intersects every conjugacy class satisfying (1) (so that r p ^ r|) and the non-split Cartan 
subgroup 



C «-={( £ °6 a)} cSi < 2 - F E 



intersects every conjugacy class satisfying (2), where e G F* is a fixed non-square element 
(so that t p ^ rj*). 

By applying Proposition 12.71 to compute the expectation and second moment on the two 
extreme sides, we find 

3 + 0{p~ l ) < E(r p ) < 3 + 0(p^), 11 + 0(p^) < E(r 2 ) < 11 + 0(p -1 )- 

which proves (1). 

To prove (2), fix some k ^ 2. We denote 

r* = max(r 1 *, r* ), r' p = max(ri, r 2 , r 3 ), 
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We have the equality of events 

{ Tp = k} = {r p = r' p = k} U {r p = k < r' p }, 
which is of course a disjoint union. Then we note that 

P(r p = k < r' p ) < Yl P ( t p = * t p > 3)- 

But clearly, if t* = j and r* < r' p , either one of the conjugacy classes (X[, . . . ,Xj) has 
trace zero, or otherwise we must have tL — T3 > j ^ 2. In the first case, since all X n have 
the same uniform distribution, the probability is at most 

j'P(Tr(X») = 0) « jp' 1 

for all p ^ 2 (again by looking at conjugacy classes for example). In the second case, we 
have 

P(r 3 > 3) < P(ts > 2) < p~ 2 . 
Combining this with the equality of events we found, it follows that for k fixed, we have 

P(r p = k) = P(r p = r p = k) + 0(p- 1 ) 

where the implied constant depends on k. 
Next we note that 

K = k } = {t p = r' p = k}U {r' k =p, r p < k}, 

again a disjoint union. As above, we find that 

fc-i 

P(r' k =p, r p < k) ^ ^P(r; = j < r' p ) « p" 1 

i=i 

where the implied constant depends on k, and hence we have finally 

P(T p = k)=P(T> = k)+0(p- 1 ), 

and the result now follows easily: first, by arguments already used, we have 

P{r' p = k) = P(max(r 1 , r 2 ) = k) + 0{p^) 

and then we are left with a coupon collector problem with two coupons of roughly equal 
probability by (15. 2p . This gives 

P(max(r 1 ,r 2 ) = k) = pt^ + P k 2 ~% = 2 (\ + 0{p~ l )f = + 0(p~ l ) 

for p ^ 2, the implied constant depending on k. □ 

Remark 5.3. Part (2) states that the waiting time tsl(2,f p ) converges in law, as p — > +00, to 
the waiting time for a coupon collector problem with two coupons of probability 1/2. Intu- 
itively, those represent finding matrices with split or irreducible characteristic polynomial. 

Remark 5.4. Recent results of Fulman and Guralnick (announced in [FG] ) should lead to a 
similar good understanding of c(G(F g )) when G is a fixed (almost simple) algebraic group 
over Q. Indeed, their results should also be applicable to situations with rank going to 
infinity, which are analogue of the symmetric and alternating groups that we consider now. 
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6. Symmetric and alternating groups 

We now come to the case of the symmetric groups & n and alternating groups A n . Here 
we have the following result, which is a precise formulation of a result essentially conjectured 
by Dixon |D1[ Abstract], following McKayQ 

Theorem 6.1. For n ^ 1, we have 

c(6 n ) x 1, c{A n ) x 1, c 2 (6 n ) x 1, c 2 (A n ) x 1. 

In fact, there exists a constant c > 1 such that, for all n ^ 1, we have 

E(c Te ") < 1, E(c TA «) < 1. 

The proof is based on the following difficult result of Luczak and Pyber, the proof of which 
involves a lot of information on symmetric groups. 

Theorem 6.2 (Luczak and Pyber). For any e > 0, there exists a constant C depending only 
on e such that 

P((X{ , ...,Xl) generate 6 n ) > 1 - e 
for all rn ^C and all n ^ 1. The same applies to A n . 

This is proved in |LPj , improving earlier work of Dixon |Dlj . 

Proof. We need only prove that the exponential moments E(c T ") are bounded for some c > 1, 
where r n = tq k with G n = & n (the A n case is similar). 
From Theorem I6.2[ there exists m ^ 1 such that 

(6.1) P({Yf, . . . , Yl) do not generate 6 n ) ^ - 

for any family of independent, uniformly distributed random variables on G n . 

Now let k ^ 1 be given; we can partition the set {1, . . . , k — 1} in \_(k — l)/mj ^0 subsets 
of size m and a remainder, and we observe that if r n = k, for each of these subsets /, we 
have 

by independence and (16.11) . Since all those sets are disjoints, we get 

/I \ l(fe-l)/ml 
P(Tn = *) < (~) < 2 l-( fc -DM 

for k ^ 1, and then, for any c ^ 1, we have 

E(c r ") = ckp ( T n = k) ^ 2 1+1/m (c2 1/m ) fc 

which converges, and is independent of n, for any c with 1 < c < 2 1 /™. □ 

In view of this, the following question seems natural: 
Question 1. Is it true or not that for all c > 1, we have 

E(c Te ") < 1 

for n ^ 1 (and similarly for A n )l 

8 This conjecture is imprecisely formulated in |D1) . where the "expected number of elements needed to 
generate & n invariably" seems to mean any r{n) for which P(c(&„) > r(n)) — > 0. 
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Another natural question, also suggested by Dixon, is: 

Question 2. Do the sequences (c((3 n )) (or (c(A n ))) converge as n — > +00? If they do, can 
their limits be computed? 

Our guess is that the answer is positive. In fact, we now present a heuristic model that 
suggests this and predicts the value of 

lim c(A n ). 

n— >+oo 

Our first step is to apply Corollary 12. 181 to a suitable "essential" set of maximal subgroups 
of symmetric groups to obtain a simpler waiting time that is asymptotically close to c(A n ) 
(or to c(& n )). The required result is again one due to Luczak and Pyber [LPj . 

Theorem 6.3 (Luczak and Pyber). For n ^ 1, let S n be the set of g e & n such that g 
is contained in a subgroup H of & n , distinct from A n , and such that G acts transitively on 
{l,...,n}. Then we have 

lim f n (S n ) = 0, 

where v n {A) = \A\/\@ n \ is the uniform density on the symmetric group. 
Corollary 6.4. For n ^ 1 and 1 ^ i < n/2, let 

Hi, n = {ge& n I g-{l,...,i} = {l,...,i}} 

be the subgroup of & n leaving {1, . . . ,i} invariant. Let H' in = H^ n D A n . Then the Hi n - 
resp. H' in - are maximal subgroups of & n - resp. A n -. Moreover, let 

M n = {A n } U {H hn I 1 < i < n/2} C max(6 n ), 
M' n = {H' in I 1 < i < n/2} C max(A n ). 



As in Proposition 2.17, let r n , resp. f' n , be the waiting time before conjugacy classes in 
each subgroup of M n , resp. M' n , has been observed. Then we have 

c(e n ) = E(f n ) + o(l), c 2 (6„) = E(f n 2 ) + o(l), 

as n — >■ +oo ; and similarly 

c{A n ) = E«) + o(l), c 2 (A n ) = E(«) 2 ) + o(l). 

Proof. It is known that the Hi^ n are (representatives of) the conjugacy classes of maximal 
intransitive subgroups of & n . Thus, we find by definition of S n that 

U ^ = s n , 

MGmax(6„)-M„ 

and hence the result follows immediately from Corollary 12.181 and Theorem 16. 3[ which pro- 
vides us with the assumption (12.191) required. □ 

In particular, in approaching Question [21 it is enough to consider the expectations and 
second moment of the random variables f n and f' n . Those are combinatorially simpler, or at 
least more explicit. 

In particular, note the following: an element a G & n is conjugate to an element of H i n C 
& n if and only if, when expressed as a product of disjoint cycles of lengths £j(c) ^ 1, 
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1 $C j ^ vj(a), say, has the property that a sum of a subset of the lengths is equal to i: for 
some J C {1, . . . , w(a)}, we have 

Note also that this applies equally to an element a in A n : the element is conjugate to 
H' in C A n if and only if the property above is true for its cycle lengths computed in & n 
(although these cycle lengths do not always characterize the conjugacy class of a in A n ). 

In particular, conjugacy classes (a{, ... , a k ) in (5^ (or A^) generate a transitive subgroup 
of & n (or A n ) if and only if n (which always occurs as the sum of all lengths) is the only 
common such sum for all o~j. (Indeed, if % < n occurs as a common subsum, we can assume 
that i ^ n/2, and then it is possible to select elements in each conjugacy class which all belong 
to H i n , so that the conjugacy classes can not generate invariably a transitive subgroup; the 
converse is also simple.) 

Now, we come to the model when n — > +00. The distribution of the set of lengths of 
random permutations is a well-studied subject in probabilistic group theory, and this allows 
us to make a guess as to the existence and value of the limit. 

Indeed, for i ^ 1, consider the map 

Wi : G n -»■ {0, 1, . . .} 

sending a to the number of cycles of length i in its decomposition as product of disjoint 
cycles (for i — 1, this is the number of fixed points; for i ^ n + 1, of course, this is zero, 
but it will be convenient for the asymptotic study to allow arbitrary i). Now consider, for 
each n ^ 1, any random variables s n , o~ n uniformly distributed on & n and A n , respectively. 
Then the following is a consequence of well-known results dating back to Goncharov |Goj : 
for fixed i, as n — > +00, the random variables Wi{a n ) converge in law to a Poisson random 
variable with parameter i.e., we have 

(6.2) lim P{wi{o- n ) = k) = e~ 1/l --J-, for fixed k ^ 0. 

Moreover, the limits for distinct values of i are independent, i.e., for any fixed finite set I 
of positive integers, we have 

lim P(vji(a n ) = ki for all i E I) = TT e -1/i — ^j-r. 

iei 1 

More precisely, this is proved (and much more precise results) for symmetric groups in, 
e.g., |AT| Th.l ] or |ABT| Th. 1.3]. The case of alternating groups can be deduced from this 
using the fact that the indicator function of A n in & n is given in terms of cycle-lengths by 

I _J_ ^ •Q-C72+'E74+— 

2 ' 

For instance, for fixed j, the characteristic function of Wk{cr n ) is 

By Goncharov's result, the first term converges for every t G R to the desired characteristic 
function of a Poisson variable with parameter 1/k; for the second term, we can use the method 
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of Lloyd and Shepp \LS\ §2]. Assuming k = 2k' is even (the other case being similar), one 
finds (see in particular |LS| (3)]) that the expectation over & n is the coefficient of z n in 

_l_ exp ^ (e ^) _ i)) Yl exp(-z 2 7j) = (1 + *)exp(y(l - e tt )), 

m 
&*? 

and since this function (of z G C) is regular at z = 1, those coefficients converge to for 
every fixed t. This computation proves (16. 2p . 

It seems therefore reasonable to use a model of Poisson variables to predict the limit of 
Chebotarev invariants of alternating groups. For this purpose, let A be the set of sequences 
of non- negative integers; we denote the i-th component of t G A by zui(£). Let v& be 
the infinite product (probability) measure on A such that the i-th component £i is distributed 
like a Poisson random variable with parameter This set A is meant to be like the set 
of conjugacy classes of an infinite symmetric group, and indeed, from the above, we see that 
for any finite I of positive integers and any fc, ^ defined for % G /, we have 

lim P{wi(a n ) = ki for all i G I) = Va{{£ G A \ tUi(£) = ki, i G /}). 

Now consider an infinite sequence (Xk)k^i of .A-valued, independent random variables, 
identically distributed according to ia We look at the following waiting time: 

r A = min{k ^ 1 | f| SiX,) = {+oo}}, 

where, for I G A, we denote by 5(f) C {0, 1, 2, . . . , } U {+00} the set of all sums 

^^ibi, where ^ hi ^ ct7j(f) 

(note the usual shift of notation from our description of the case of fixed n: the sequence of 
lengths of cycles occuring in a permutation is replaced by the sequence of multiplicities of 
each possible length). Then our guess for the limit of c(A n ) is that 

lim c(A n ) = E(ta). 

n— >+oo 

We hope to come back to this question in a future work. 

7. NON-ABELIAN GROUPS: NUMERICAL EXPERIMENTS 

In this section, we give some tables of values of the Chebotarev invariant (and the sec- 
ondary invariant) for some non-abelian finite groups. Although those are clearly rational 
numbers, we list real approximations only because the "height" of those rationals grows very 
fast. 

The computations are feasible even for fairly large and complicated non-abelian groups, 
because they may have few conjugacy classes of maximal subgroups, and not too many 
conjugacy classes. For instance, the Weyl group W(Es) (one of our motivating examples) 
has 9 conjugacy classes of maximal subgroups, and 112 conjugacy classes. If these data 
are available to suitable software packages, Proposition 12.71 provides a way to compute the 
Chebotarev invariants, though this is at best an exponential-time algorithm (due to the 
necessity to sum over all subsets of max(G)). 
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The computations here were done for the most part with Magma (see [M]), using the script 
included in the Appendix. The correctness of the results was checked partly by independent 
computations with the open-source package GAP (see |GAPj ). and by checking that the 
results agree, for cyclic groups and groups F p , with the theoretical formulas of the Sectional 
They are also in good agreement, in the case of PSL(2,F P ), with the asymptotic result of 
Section [5j Hence, altogether, we have very high confidence in these values. 

The computations were relatively fast; usually there was a sharp threshold between com- 
puting for one group in a family in less than an hour, and the next one proving infeasible due 
to the exponential growth of the number of subsets of max(G). As an indication of timing, 
the computation for PSL(Q, F3) with Magma (version 2.14.15) took about 42 seconds on a 
2.5 GHz Core 2 processor. 

Below, we include tables for the alternating groups A n , for the symmetric groups & n , for 
the groups PSL(2, F p ) with p prime =Sj 150 (though the computations can be done for p 
quite a bit larger, we do not include the results which are not particularly enlightening), for 
PSL(3,F P ), PSL(A,F P ), PSL(n, F 2 ), PSL(n,F 3 ), PSL(n, F 4 ), PSL(2,F 2 »), P£L(3,F 2 »), 
Sp(2g, F3). (Note that, in general, the computations tend to run quite a bit faster for 
simple groups.) We also include a table of the "partial" invariants E(f^) and E((f^) 2 ) of 
alternating groups defined in Corollary 16.41 Note that although we have shown that these 
are asymptotically converging to the Chebotarev invariants themselves, the convergence is 
by no means visible! There is also a table for the Borel subgroup of SL(3, F p ), namely 

(x r s\ 
y t) I (r,s,t)eF 3 p , (x, y, z) E (F*) 3 , xyz = l}, 
zj } 

as another example of a solvable group. 

Another table lists some more "sporadic" groups; the names of those groups in the table 
should be self-explanatory. For instance, D 2n is the dihedral group of order 2n, W(R) denotes 
the Weyl group of a root system of type R; Sz(8) and Sz(32) are Suzuki groups. Sporadic 
simple groups are named in a standard way: 

• Mathieu groups: M n , where n E {11, 12, 22, 23, 24}; 

• Janko groups: Jk, where k E {1,2,3,4}; 

• Second and thirs Conway groups: C02, C03 (the first Conway group is too big); 

• Tits group T; 

• MacLaughlin group McL; 

• Higman-Sims group HS; 

• Helde group He. 

The group Rub at the end of the table is the Rubik's group (the subgroup of ©43 that gives 
the possible moves on the Rubik's Cube; computing c(Rub) takes about two days on a fast 
Opteron; this group has 20 conjugacy classes of maximal subgroups and 81120 conjugacy 
classes). In order to ease checking, the url 

http : / /www . math . ethz . ch/~kowalski/other-groups . mgm 

contains a Magma file where each group in this list is constructed explicitly. 

It also possible to exploit the databases of small groups, or of transitive groups, or primitive 
groups, to compute the Chebotarev invariants for, say, all groups of a given small order (up 
to isomorphism), or for all transitive permutations groups of small degree. The latter is of 
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course particularly interesting from the point of view of Galois theory, and the groups F q x F* 
which appear as transitive permutation groups of degree q (and in Galois theory as Galois 
groups of Kummer extensions of prime-power degree, i.e., splitting fields of polynomials of 
the type X q — a) are very noticeable, having much higher Chebotarev invariants than the 
other groups despite their rather small order (see the example in the table for transitive 
groups of degree 17 - noting that the group with Chebotarev invariant roughly 8.88 is the 
index 2 subgroup of Hn denoted Ci in Section H|). We include a figure of the empirical 
distribution of values for the Chebotarev waiting time for H 3 i (chosen because q — 1 = 30 
has "many" divisors). 

We also include a figure with an histogram showing the distribution of the Chebotarev 
invariant for the 840 distinct groups of order 720 (up to isomorphism). Note that this data 
also indicates that the invariant is far from injective (as can be guessed from its dependency 
on relatively little data): there are only 188 distinct values of c(G) for \G\ = 720; the value 

469589438194474533813031879 

~ 5.836158 .... 

80462083849550829871525080 

occurs with maximal multiplicity (it arises 39 times). 

Note that for simple groups (or groups which are nearly so), the relation between c(G) 

and 02(G) seems relatively regular, but there is certainly no strict monotony in terms of the 

order; see, e.g., the cases of alternating groups A n , where sorting according to the value of 

c(A n ) leads to the following rather bizarre ordering of the segment 2 «C n ^ 21: 

2, 3, 13, 19, 17, 11, 5, 10, 14, 20, 21, 16, 18, 15, 4, 6, 12, 9, 7, 8 ; 

the ordering with respect to C2(A n ) is slightly different, namely: 

2, 3, 13, 19, 11, 17, 10, 14, 21, 20, 16, 18, 15, 5, 12, 6, 9, 4, 7, 8. 

And the orderings for c((3„) and C2(@ n ) are also different: 

2, 3, 7, 11, 13, 9, 17, 19, 5, 15, 21, 16, 20, 4, 14, 18, 12, 10, 8, 6, 

and 

2, 3, 7, 11, 13, 9, 17, 19, 15, 5, 21, 16, 20, 14, 18, 12, 8, 10, 4, 6, 

respectively. Note however that in Table 2, if we fix the parity of n, the invariants ~Ei(f! 2n ) and 
E(f"2„ +1 ) seem monotonically increasing. This indicates that they are indeed very natural 
objects to study. 

8. Arithmetic considerations 

In this short section, we indicate the (expected) number-theoretic connections of our work. 

First, let K be a Galois extension of Q with group G. For each prime p that is unramified 
in K, we have a well-defined Frobenius conjugacy class Fr Pi x £ G*. For simplicity, we denote 
F r p,if = 1 when p is ramified in K. The Chebotarev density theorem says that 

(8.1) Ihn lfr<» = ?y = C}l = M 

where C G G^ is a fixed conjugacy class of G and 7r(y) is the usual prime-counting function, 
i.e., the number of primes p ^ y. 

Now fix a real number y large enough, so that every conjugacy class of G is of the form 
Fr Pt K f° r some p ^ y. For each i ^ 1, select uniformly and independently a random prime p 
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Figure 1 . Empirical distribution of the waiting time for H^i 
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Figure 2. Distribution of the Chebotarev invariant for groups of order 720 
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from the set {p : p ^ y} and define X\ y = Fr Pj ^- We thus have a sequence of independent 

and identically distributed random variables X(y) = (Xf ) in GK As usual, we define the 
waiting time 

Tx( y ),G = min{n ^ 1 | (X[ >y , . . .,Xl >y ) generate G} G [1, +00]. 
Using the Chebotarev density theorem, one obtains easily 

lim E(rx(y), G ) = c(G). 

Therefore, in an imprecise way, c(G) can also be thought of as the expected number 
of "random" primes p needed for Fr Pi ^ to generate G = Gal(/i/Q). Indeed, this is our 
motivation for using the name "Chebotarev invariant" . 

Of course in practice, one usually considers the (non-random) sequence Fr 2j x, Fr 3i ^, Fr 5! ^, 
Fr? t K, ■ ■ ■■ We now explain, informally, what can be expected to happen in that situation. 
The deterministic analogue of the Chebotarev waiting time is given 

t(K) = mm{k ^ 1 | the first k conjugacy classes Fr 2i x5 ■ • • , ^ T p k ,K generate G}, 

where p^ is the k-th prime number. 

However, for a fixed K/Q, the value of t(K) might diverge considerably from c(G). So 
we suppose we have some family % of finite Galois extensions of Q (or another base field), 
all (or almost all) of which have Galois group Gal(-fT/Q) ~ G, a fixed finite group, and that, 
for all values of some parameter x ^ 1, we have finite subfamilies % x (which exhaust % 
as x — > +00) and some averaging process for invariants of the fields in X, denoted E x (for 
instance, one might take 

E.(«(iO) = mH E 

' Kex x 

but other weights, involving multiplicities, etc, might be better adapted). 

Using this, we can define Chebotarev invariants for the family % by averaging: 

c(X x ) = E x (r(K)), c 2 (X x ) = E x (t(K) 2 ). 

The basic arithmetic question is then: for a given family, is it true that c(X x ) is, for x suf- 
ficiently large at least, close to c(G) (and similarly for the secondary Chebotarev invariant)? 
The basic reason one can expect this is the Chebotarev density theorem (18.11) . We want to 
point out a few difficulties that definitely arise in trying to make this precise. 

First of all, quantifying the Chebotarev density theorem is hard: it almost immediately 
runs into issues related to the Generalized Riemann Hypothesis; even in the seemingly trivial 
case where G = Z/2Z (quadratic extensions of Q), the basic question of estimating the size 
of the smallest prime p for which the corresponding Frobenius is non-trivial, i.e., the smallest 
quadratic non-residue modulo p, in terms of the discriminant of the field, is unsolved (see, 
e.g., Prop. 5.22, Th. 7.16] for conditional and uncontional results in that case). This is 
a problem because if we sum with uniform weight, a single "bad" field Kq can destroy any 
chance of approaching the Chebotarev invariant. Indeed: note that in that case 

(8.2) E x (t{K)) > -±-k mm (K ) 

where 

kmin(K) = mm{k ^ 1 | Fr Pi ^ ^ 1} 
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is the index of the first non-trivial Frobenius conjugacy class. In the current state of knowl- 
edge, it can be that there exists Kq with 

(K ) > disc(K ) A 

for some constant A > (see |LMOj ); on the other hand, if the family % is defined as 
that of splitting fields of monic polynomials of degree n, and the subfamily % x is that of 
polynomials of height ^ x, then we know that most K G % have Galois group & n , that 
\% x \ = (2x + l) n if x is an integer, and the discriminant is obviously often also at least a 
power of x. Thus (I8.2p might already be bad enough to preclude any comparison. On the 
other hand, on the Riemann Hypothesis, we have 

(K) < (log disc (TT)) 2 , 

(where the implied constant depends on G), and the problem would then be alleviated. 
Another issue is that one can not expect, as stated, to have 

lim c(X x ) = c(G) 

for interesting families for the simple reason that the statistic of small primes is typically 
not the uniform one, i.e., if we fix a prime p, we can not expect to have 

lim E X (1 {F =cS} ) = z/ G (c a ), 

even if we assume that all the fields involved are unramified at p. 
For instance, consider % the set of cubic polynomials 

X 3 + a 2 X 2 + a x X + oq, 

with di G Z, with % x those where |aj| ^ x for % = 0, 1, 2, and count them uniformly. Take 
p = 5 and consider only polynomials with no repeated root modulo 5 and splitting field of 
degree 6; then, asymptotically, the conjugacy Frobenius at 5 will be distributed in ©3 = G 
as dictated by the factorization of the polynomial modulo 5. One finds easily that there 
are 100 monic cubic polynomials in F5LY] with non-zero discriminant (there are 25 with 
repeated roots), among which: 

• 10 split in linear factors, i.e., a density 1/10; 

• 40 are irreducible, i.e., a density 4/10; 

• 50 split as a product of one linear factor and one irreducible quadratic factor, i.e., a 
density 1/2. 

This is in sharp constrast with the density of the three corresponding conjugacy classes in 
63, which are respectively: 

• 1 /6 for the identity class; 

• 1/2 for the 3-cycles; 

• 1/3 for the transposition. 

In particular, not even the relative frequencies are preserved! On the other hand, it is well- 
known that if p is increasing, the discrepancy between the distribution of the factorization 
patterns of squarefree polynomials modulo p and the density of conjugacy classes disappears: 
we have 

— 1{/ G F P [X] I / squarefree of degree n with Frj = c"}| ~ ^g(c') 
uniformly for all conjugacy classes c" G G = & n . 
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This suggests that it is likely that one can prove some relevant results: one would consider 
some increasing starting point s(x) ^ 2 and a modified waiting time 

t x (K) = min{A; | the first k conjugacy classes Fr Pi ^ with p ^ s(x) generate G} 

and hope to prove (possibly under the Generalized Riemann Hypothesis, possibly uncondi- 
tionally after throwing away a few "bad" fields) that 

lim E x (t x (K)) = c(G), 

x— >+oo 

for suitable s(x). One may guess that for polynomials of height ^ x and fixed degree n (and 
G = & n ), this would be true with s(x) x logo;. 

9. Remarks and problems 

We finish with a few more remarks and problems. 

• (What does the invariant "know"?) As a bare numerical invariant of a finite group, 
the Chebotarev invariant seems to be fairly subtle. For instance, we see from Section[3] 
that it "knows" that vector spaces over finite fields are in some sense very similar 
for varying base field, but that they become also "simpler" as the cardinality of the 
base field grows. It also seems to know that non-reductive finite matrix groups are 
worse-behaved than reductive ones (as shown by the results for H q ). What else does 
the invariant reveal? 

• (A method for upper bounds) There are, in the literature, quite a few results about a 
finite group G of the type: "if a subgroup H contains elements in some set C\, some 
other set C2, some other set C m , of conjugacy classes, then H is in fact equal 
to G n . For instance, a lemma of Baer quoted by Gallagher [G] Lemma, p. 98] says 
that there is no proper subgroup H of & n which (1) contains an n-cycle, (2) contains 
a product of a transposition and cycles of odd lengths, (3) contains an element of 
order divisible by a prime p > n/2. Another such result is the Lemma [5.21 of Serre 
for SL(2,¥ p ), and we also mention [JKZl Lemma 3.2] for another example with the 
Weyl group W(E 8 ), and there are many other such results used, e.g., for proving 
concrete cases of Hilbert's Irreducibility Theorem. 

With this notation, and assuming we work with a sequence of independent and 
uniformly distributed G- valued random variables (X n ), this means that we have 

ra «S T Cl ,...,c m = max(r Ci , l^j^m), 

where 

r Cj = mm{n ^ 1 | XJgQ}. 
From Proposition 12.91 we obtain easily an upper bound 

(-1 )\ T \ +1 

(9-1) c(G) ^ E(r C: <;. ) = 



0^JC{l,...,m} ^(U fe/Cj 



and one may hope to approximate c(G) by choosing wisely the sets (Cj). 

However, it is not clear at all to what extent this can approach the truth. Here are 
some examples: 

(1) Baer's lemma gives only an upper bound 

c(6 n ) < n 
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as n — > +00, which is quite weak compared with Theorem 16.11 (it is dominated by 
the time required to obtain an n-cycle). How far is this from the best possible result 
that can be obtained in this way, and how far is the latter from Theorem 16. 11 / 
(2) Consider G = F 2 with p odd. It is possible to take 

Ci = {(x, y) G F 2 — {0} I y 7^ and xy~ l is a square in F p }, 
C 2 = F 2 p -{0}-C 1 . 

The point is that whenever (v , w) G C\ x C 2 , w and v are not on the same line 
through the origin, so (v, w) generate F 2 . Since 

|Ci| = |C 2 | = (p 2 - l)/2, |d U C 2 \ = P 2 -1, 



this leads to 



/t,^ 11 1 3» 2 



HCi) v(C 2 ) v{Cx U C 2 ) p 2 -l' 
which asymptotically requires one more step on average than the right Chebotarev 
invariant (given by (E2])), namely c(F 2 p ) = (2p 2 + p)/(p 2 - 1). It seems also that this 
type of sets is essentially best possible for applying this upper bound in this case. 

(3) Consider G = W(Es), the Weyl group of Eg. There is a non-trivial homomor- 
phism 

e : W(E 8 ) -> {±1}, 

and in |JKZ| Lemma 3.2], jointly with F. Jouve, we proved that one could take 
Ci = ker(e), C 2 the union of the conjugacy classes of w and w 2 , where w G W(E 8 ) is 
a Coxeter element; the density of C 2 is 1/15 and we then get the upper-bound 

2 + 15 - 30/17 = 25.23... 

instead of the correct value 4.194248 .... 

(4) For SL(2, F p ), Theorem 15 . 1 1 shows that the sets C\, C 2 , C3 given by Lemma [531 
give an asymptotically optimal answer (and this is an essential ingredient in the 
proof). 

Despite this relative inefficiency, it is interesting to notice that in applications of 
sieve methods to probabilistic Galois theory (as was the case in P) and [JKZ]0 it is 
this type of distinguishing families which can be used in estimating how rare "small" 
Galois groups are in certain families, and in fact it is the quantity 

m 

which occurs naturally as coefficient in a "saving factor" of the large sieve; see, 
e.g, |K1| p. 57], where the question of minimizing this by varying the sets was raised 
explicitly for symmetric groups. 

(General upper bounds?) A first problem is to bound c(G) from above, in a mean- 
ingful way. Since we have 

Memax(G) 



only implicitly in the latter. 
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we obtain 



MGmax(G) V ' 

from ( 12. 3ft . Together with ( 12.15p . this gives an upper bound 

c(G)<|G| E M 

Memax(G) ' ' 

which is close to being sharp for the groups H q of Section HJ indeed, if q is odd, then 
Lemma 14.21 gives 

Kemax(ff 9 ) 1 1 ^ i\q-l K ' 

£|qr — 1 2<£|q-l 

(where £ runs over prime divisors of q — 1). \l q = 2i + 1 (£ odd prime) is a Sophie 
Germain prime, this gives 

9-1 3(g + l) 



1^1 



F A 

Kemax(/f,) 



g+2+ E ^ = 9 + 

2<%-l 

which is off, asymptotically, only by a factor 3/2 from the value 

c{H q ) = q + 0{q~ l ) 

that follows from ( 14. ip . Of course, it is not known that there are infinitely many 
Sophie Germain primes, but for q = llyli + 1, with li prime, we have 

J_ |g + £i+^ 2 + 2 Mtx + h 
|^|~\g + £i + 2 if ^ =4 

SC 2g. 

By sieve methods, it is known that there are infinitely many primes q for which 
either q is a Sophie Germain prime, or is 2£i£ 2 + 1, and hence one sees that the 
"trivial" estimate ( 19. 3 p above can not be improved by more than a constant in full 
generality. On the other hand, it is very far off for many groups: for a random 
example, it gives 

4.7820 . . . = c(A 7 ) ^ 93. 

It would be more interesting to have a decent upper bound in terms of the order 
of G only. Here, using the set of all conjugacy classes in (19. ip . we get as an upper 
bound from the contribution of singletons that 

(where Nc{g) is the normalizer of g in G). This gives trivially 

c(G) ^ |Gf , 
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but this seems unlikely to be close to the truth (for 1). For instance, since 

c(H q ) = q ~ y/\Hj, 

one may wonder if H q is also (essentially) extremal in this sense, i.e., one may ask 
whether an estimate 

c(G) < ^/\G\ 

holds for all G. (Certainly for |G| = q(q — 1) with q ^ 43 prime, it is experimentally 
true that H q maximizes the Chebotarev invariant). 
• (Other classes of groups?) There are many classes of groups for which it should be 
possible to understand the behavior of the Chebotarev invariant, at least asymptoti- 
cally. For instance, one can consider non-reductive subgroups of finite matrix groups, 
e.g., the standard Borel subgroup (upper triangular matrices) of GL(n,F q ). In fact, 
solvable groups seem particularly interesting. 

Appendix: Magma script 

The following script can be used to compute the Chebotarev invariant (and the secondary 
invariant) using Magma, by applying the formulas (I2.3P and (12. 4p . The output is given as 
real approximations since usually the denominators are unwieldy. Also note that because 
of the use of the construct Subsets (M), this script only applies to groups with at most 29 
conjugacy classes of maximal subgroupsj^ to - try to - compute further, one would have to 
replace the loop over subsets obtained in this manner with a hand-rolled one. 

A similar GAP script is available upon request, as well as a Sage version, which basically 
calls the GAP group theory routines. However, these versions are much slower. 

The last routine in the script is useful for "empirical" study of the probabilistic model. 

// The following calculates J such that 

//J [k] [i] =true if and only if the k-th maximal subgroup 
// of G intersects the i-th conjugacy class of G 

MCIntersectionMatrix:=f unction(G,C,f ,M) 

J : = [ [false : i in [l..#C]] : k in [1 . . #M] ] ; 
for k in [1. .#M] do 
H := M[k] 'subgroup; 
CH := ConjugacyClasses(H) ; 
for j in [1. .#CH] do 

J[k] [f (CH[j] [3])] := true; 
end for; 
end for; 
return J; 
end function; 

// This returns [c,s] where c is the Chebotarev invariant of G 
// and s the secondary invariant. 

Chebotarev: = function (G) 
if IsTrivial(G) then 

10 For alternating groups A ni this means n ^ 47, or n <E {49, 51, 53}. 

36 



return <1.0,1.0>; 
end if ; 

C := ConjugacyClasses(G) ; 

f := ClassMap(G) ; 

M := MaximalSubgroups(G) ; 

J := MCIntersectionMatrix(G,C,f ,M) ; 

c:=0.0; s:=0.0; 

for I in Subsets ({1 . .#M}) do 
if #1 ne then 
v:=0; 

for i in [1. .#C] do 

if forall(t) {k: k in I I J[k][i]} then 
v:= v + C[i] [2]/#G; 

end if ; 
end for; 

c := c + (-l)~(#I+l)/(l-v) ; 

s := s+ (-1)~ (#I)/(l-v)*(l-2/(l-v)) ; 
end if ; 
end for; 
return ( [c , s] ) ; 

end function; 

// Compute empirical Chebotarev invariant. 
// The optional parameter steps is the number 
// of iterations to do. Example: 
// > EmpiricalChebotarev(Alt (7) : steps: =10000) ; 

EmpiricalChebotarev: =f unction(G : steps :=1) 
total :=0; 

C := ConjugacyClasses(G) ; 
f :=ClassMap(G) ; 
M := MaximalSubgroups(G) ; 
J := MCIntersectionMatrix(G,C,f ,M) ; 
for count in [1.. steps] do 
nb : =0 ; 

vprint Userl: "Iteration, ", count; 
// Start with all subgroups 
possible:=[ 1. .#M ] ; 
while possible ne [] do 

g : =Random(G) ; 

nb := nb+1; 

index :=f (g) ; 

// Only those subgroups containing the class of g remain 
possible: =[ k : k in possible I J [k] [index] ]; 

end while; 

total : =total+nb ; 
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end for; 

return total/steps, total/steps*l . 0; 
end function; 
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Table 1. Chebotarev invariants of A. 



n 


Order 


/ A \ 

c(A n ) 


/ A \ 


2 


1 


1.000000. . . 


1.000000. . . 


3 


3 


1.500000. . . 


3.000000. . . 


4 


12 


4.409091. . . 


29.71074. . . 


5 


60 


4.136364. . . 


22.64463. . . 


6 


360 


4.439574. . . 


25.49003. . . 


7 


2520 


4.782001. . . 


29.98671. . . 


8 


20160 


4.939097. . . 


31.98434. . . 


9 


181440 


4.637463. . . 


26.35009. . . 


10 


1814400 


4.145282. . . 


21.73709. . . 


11 


19958400 


4.092974. . . 


21.08692. . . 


12 


239500800 


4.444074. . . 


24.14188. . . 


13 


3113510400 


4.016324. . . 


20.51475. . . 


14 


43589145600 


4.212753. . . 


22.16514. . . 


15 


653837184000 


4.289698. . . 


22.51291. . . 


16 


10461394944000 


4.239141. . . 


22.21416. . . 


17 


177843714048000 


4.089704. . . 


21.12890. . . 


1 8 
lo 


OZVl lo0o0zo04UUU 


A 9/1 81 QQ 


//.ooUoo. . . 


19 


60822550204416000 


4.072274. . . 


21.08656... 


20 


1216451004088320000 


4.229094. . . 


22.20516... 


21 


25545471085854720000 


4.238026. . . 


22.19523... 


22 


562000363888803840000 


4.240513... 


22.33370. . . 


23 


12926008369442488320000 


4.131077... 


21.54514... 


24 


310224200866619719680000 


4.282667. . . 


22.58460. . . 
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Table 2. "Partial" Chebotarev invariants of A 



n 


Order 


E«) 


E(«) 2 ) 


3 


3 


1.500000. . . 


3.000000. . . 


4 


12 


2.123377. . . 


5.874009. . . 


5 


60 


2.500000. . . 


10.00000. . . 


6 


360 


2.649424. . . 


9.187574. . . 


7 


2520 


3.243247. . . 


16.47701. . . 


8 


20160 


2.812743. . . 


10.71084. . . 


9 


181440 


3.133704. . . 


13.97383. . . 


10 


1814400 


3.115450. . . 


13.08967. . . 


11 


19958400 


3.399573. . . 


15.88920. . . 


12 


239500800 


3.225496. . . 


14.16483. . . 


13 


3113510400 


3.402011. . . 


15.56383. . . 


1 A 


4oOoy 1400UU 


O.OO ( DDI. . . 


10. lO 1 4Z. . . 


15 


653837184000 


3.504050. . . 


16.37350. . . 


16 


10461394944000 


3.385358. . . 


15.32752... 


17 


177843714048000 


3.544719. . . 


16.55867... 


18 


3201186852864000 


3.497980. . . 


16.21775... 


19 


60822550204416000 


3.625919... 


17.22183... 


20 


1216451004088320000 


3.530703. . . 


16.46076. . . 
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Table 3. Chebotarev invariants of © 



n 


Order 


c(©„) 


C 2 (©n) 


2 


2 


2.000000. . . 


6.000000. . . 


3 


6 


3.800000. . . 


19.32000. . . 


4 


24 


4.498380. . . 


25.91538. . . 


5 


120 


4.331526. . . 


23.50351. . . 


6 


720 


5.610738. . . 


37.63260. . . 


7 


5040 


4.115230. . . 


21.20184. . . 


8 


40320 


4.626289. . . 


25.71722. . . 


9 


362880 


4.250355. . . 


22.49197. . . 


10 


3628800 


4.624666. . . 


25.76898. . . 


11 


39916800 


4.173683. . . 


21.86294. . . 


12 


479001600 


4.583705. . . 


25.11338. . . 


13 


6227020800 


4.213748. . . 


22.21319. . . 


14 


87178291200 


4.508042. . . 


24.57963. . . 


15 


1307674368000 


4.365718. . . 


23.39257. . . 


16 


20922789888000 


4.461633. . . 


24.12713. . . 


17 


355687428096000 


4.282141. . . 


22.79488. . . 


1 8 
lo 


04UZO lot UO ( ZoUUU 


A t^QI 78/1 
4.001 I o4. . . 


OA R7fiSD 


19 


121645100408832000 


4.308469. . . 


23.01145... 


20 


2432902008176640000 


4.497047. . . 


24.37207. . . 


21 


51090942171709440000 


4.391209... 


23.61488... 


22 


1124000727777607680000 


4.477492. . . 


24.29632. . . 


23 


25852016738884976640000 


4.352364. . . 


23.37533. . . 


24 


620448401733239439360000 


4.523388. . . 


24.57409. . . 



Table 4. Chebotarev invariants of transitive groups of degree 17 



Name 


Order 


c(G) 


c 2 (G) 


Z/17Z 


17 


1.062500... 


1.195312... 


C 8 C Hn 


34 


3.094697. . . 


11.81350... 


C4 C Hyj 


68 


4.890000. . . 


35.53580. . . 


Ci C Hyj 


136 


8.880953. . . 


138.3764. . . 




272 


17.21053... 


562.3851... 


PSL(2,F 16 ) 


4080 


3.200912. . . 


12.73727. . . 


7 


8160 


4.055261... 


20.84364. . . 


8 


16320 


4.067118... 


20.58582... 


A l7 


177843714048000 


4.089704. . . 


21.12890... 


©17 


355687428096000 


4.282141... 


22.79488. . . 
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Table 5. Chebotarev invariants of PSL(3,F P ) 



p 


Order 


c(PSL(3,F p )) 


c 2 (PSL(3,F p )) 


2 


168 


4.653153... 


29.48762. . . 


3 


5616 


3.845890. . . 


20.67132. . . 


5 


372000 


3.629464. . . 


18.36114. . . 


7 


1876896 


3.661481... 


18.91957. . . 


11 


212427600 


3.527819... 


17.29354. . . 


13 


270178272 


3.546344. . . 


17.55063. . . 


17 


6950204928 


3.511708. . . 


17.12456. . . 


19 


5644682640 


3.521753... 


17.25893. . . 


23 


78156525216 


3.506462. . . 


17.06878. . . 


29 


499631102880 


3.504076. . . 


17.04348. . . 


31 


283991644800 


3.508213... 


17.09800. . . 


37 


1169948144736 


3.505795. . . 


17.06906. . . 


41 


7980059337600 


3.502051... 


17.02191. . . 



Table 6. Chebotarev invariants of PSL(4, F p ) 



p 


Order 


c(PSL(A,F p )) 


c 2 (PSL(A,F p )) 


2 


20160 


4.939097. . . 


31.98434... 


3 


6065280 


4.191257... 


23.35082. . . 


5 


7254000000 


3.768197... 


18.89633... 


7 


2317591180800 


3.613602... 


17.31973... 


11 


2069665112592000 


3.530797. . . 


16.44109... 


13 


12714519233969280 


3.513963... 


16.24990... 



Table 7. Chebotarev invariants of PSL(n, F 2 ) 



n 


Order 


c(PSL(n,F 2 )) 


c 2 (PSL(n,F 2 )) 


2 


6 


3.800000. . . 


19.32000... 


3 


168 


4.653153... 


29.48762. . . 


1 


20160 


4.939097. . . 


31.98434... 


5 


9999360 


4.238182... 


25.64374. . . 


6 


20158709760 


4.456089. . . 


27.20052. . . 


7 


163849992929280 


4.335957. . . 


26.54874. . . 


8 


5348063769211699200 


4.465723. . . 


27.53266. . . 


9 


699612310033197642547200 


4.460433. . . 


27.64706. . . 
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Table 8. Chebotarev invariants of PSL(n,F 3 ) 



n 


Order 


c{PSL(n, F 3 )) 


C2{PSL{n,F 3 )) 


2 


12 


4.409091. . . 


29.71074. . . 


3 


5616 


3.845890. . . 


20.67132... 


4 


6065280 


4.191257... 


23.35082. . . 


5 


237783237120 


3.949889. . . 


21.81110... 


6 


21032402889738240 


4.123378... 


23.06449. . . 


7 


67034222101339041669120 


4.066340. . . 


22.81370... 



Table 9. Chebotarev invariants of PSL(n, F 4 ) 



n 


Order 


c{PSL{n,F 4 )) 


c 2 (PSL(n,F 4 )) 


2 


60 


4.136364. . . 


22.64463. . . 


3 


20160 


4.399979. . . 


26.39681... 


4 


987033600 


3.770618... 


19.19928... 


5 


258492255436800 


3.838194... 


20.33428. . . 


6 


361310134959341568000 


4.002927. . . 


21.57223... 



Table 10. Chebotarev invariants of PSL (2, F 2 n) 



2 n 


Order 


c(PSL(2,F 2 n)) 


c 2 (PSL(2,F 


2 


6 


3.800000. . . 


19.32000. . 


4 


60 


4.136364. . . 


22.64463. . 


8 


504 


3.437879. . . 


14.95188.. 


16 


4080 


3.200912... 


12.73727. . 


32 


32736 


3.096876. . . 


11.82191. . 


64 


262080 


3.048732. . . 


11.40623. . 


128 


2097024 


3.023623. . . 


11.19773.. 


256 


16776960 


3.011765... 


11.09826.. 


512 


134217216 


3.005965. . . 


11.04945. . 



2")) 



Table 11. Chebotarev invariants of PSL (3, F 2 «) 



2 n 


Order 


c(PSL(3,F 2 n)) 


c 2 (P,SL(3,F 2 n)) 


2 


168 


4.653153... 


29.48762. . . 


4 


20160 


4.399979. . . 


26.39681... 


8 


16482816 


3.551417... 


17.54363. . . 


16 


1425715200 


3.549690. . . 


17.47208. . . 


32 


1098404364288 


3.503357... 


17.03581... 
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Table 12. Chebotarev invariants of the Borel subgroup of SL(3, F p ) 



V 


Order 


c(B 3 (F p )) 


C2(B 3 (F P )) 


2 


8 


3.333333. . . 


13.55556... 


3 


108 


5.074442. . . 


31.76009... 


5 


2000 


7.686557. . . 


81.14365... 


7 


12348 


10.07528... 


150.8724. . . 


11 


133100 


16.38777. . . 


402.7223. . . 


13 


316368 


18.85106... 


551.0363... 


17 


1257728 


25.31072... 


978.0196... 


19 


2222316 


27.79352. . . 


1204.483. . . 


23 


5888828 


34.28491. . . 


1805.763... 


29 


19120976 


43.27249. . . 


2885.634. . . 


31 


26811900 


45.75644. . . 


3268.081... 


37 


65646288 


54.75057. . . 


4678.007. . . 


41 


110273600 


61.26132... 


5801.515... 


43 


140250348 


63.74680. . . 


6339.956. . . 
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Table 13. Chebotarev invariants of PSL(2,F P ), p ^ 150 



p 


Order 


c(PSL(2,F p )) 


c 2 (PSL(2,F p )) 


2 


6 


3.800000. . . 


19.32000... 


3 


12 


4.409091... 


29.71074... 


5 


60 


4.136364. . . 


22.64463. . . 


7 


168 


4.653153... 


29.48762. . . 


11 


660 


3.981397. . . 


20.76193... 


13 


1092 


3.293965... 


13.63659... 


17 


2448 


3.264353. . . 


13.20732... 


19 


3420 


3.259202... 


13.08533... 


23 


6072 


3.136600... 


12.18536... 


29 


12180 


3.115633... 


11.99619... 


31 


14880 


3.111661... 


11.92578... 


37 


25308 


3.088522... 


11.75723... 


41 


34440 


3.098342. . . 


11.78358... 


43 


39732 


3.071689... 


11.61064... 


47 


51888 


3.065454. . . 


11.55651... 


53 


74412 


3.060208. . . 


11.51103... 


59 


102660 


3.051900... 


11.43952... 


61 


113460 


3.051897... 


11.43943... 


67 


150348 


3.045600. . . 


11.38545... 


71 


178920 


3.046777. . . 


11.38343... 


73 


194472 


3.042989. . . 


11.36306... 


79 


246480 


3.043013... 


11.34889... 


83 


285852 


3.036689. . . 


11.30930... 


89 


352440 


3.036100... 


11.30056... 


97 


456288 


3.031998... 


11.26935... 


101 


515100 


3.032463. . . 


11.26755... 


103 


546312 


3.030308. . . 


11.25228... 


107 


612468 


3.028370. . . 


11.23855... 


109 


647460 


3.029877. . . 


11.24644. . . 


113 


721392 


3.028016... 


11.23330... 


127 


1024128 


3.024393. . . 


11.20309... 


131 


1123980 


3.024148. . . 


11.19945... 


137 


1285608 


3.022889. . . 


11.19063... 


139 


1342740 


3.022686. . . 


11.18747... 


149 


1653900 


3.020586. . . 


11.17269... 
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Table 14. Chebotarev invariants of some other groups 



Name 


Order 


c(G) 


c 2 (G) 


W(G 2 ) = D 12 


12 


4.31515...= 717/165 


23.45407. . . 


W(C A ) 


384 


4.864890. . . 


29.10488. . . 


W{F A ) 


1152 


5.417656... 


35.12470... 


GL(2,F 7 ) 


2016 


3.767768. . . 


17.29394. . . 


A 5 x A 5 


3600 


5.374156. . . 


35.41628. . . 


W(C 5 ) 


3840 


4.863533. . . 


28.13517. . . 


Mn 


7920 


4.850698. . . 


29.72918. . . 


GL(3,F 3 ) 


11232 


4.110394. . . 


22.77077. . . 


G 2 (F 2 ) 


12096 


5.246204. . . 


34.24515. . . 


Sz(8) 


29120 


3.101639... 


11.92233. . . 


W(C 6 ) 


46080 


5.792117... 


39.56093. . . 


W(E 6 ) 


51840 


4.470824. . . 


23.93050. . . 


Sp(4, F 3 ) 


51840 


4.401859. . . 


24.03143. . . 


PGL(3,F 4 ) 


60480 


3.763384. . . 


19.49865. . . 


M 12 


95040 


4.953188. . . 


29.53947. . . 


Ji 


175560 


3.423739. . . 


14.76364. . . 


M 22 


443520 


4.164445. . . 


22.70981. . . 


■h 


604800 


4.031298. . . 


19.07590. . . 


W(C r ) 


645120 


4.632612... 


25.54504. . . 


PSp(6, F 2 ) 


1451520 


5.270439. . . 


34.84139. . . 


W(E 7 ) 


2903040 


5.398250... 


36.04850. . . 


G 2 (F 3 ) 


4245696 


4.511630... 


24.06106. . . 


M 23 


10200960 


4.030011. . . 


20.98580. . . 


W(C 8 ) 


10321920 


4.928996. . . 


28.53067. . . 


T 


17971200 


4.963701. . . 


32.54160. . . 


Sz(32) 


32537600 


2.755449. . . 


9.107751. . . 


HS 


44352000 


4.484432. . . 


25.68549. . . 


■h 


50232960 


4.304616. . . 


23.42082. . . 


W{Cg) 


185794560 


4.716359... 


26.41344. . . 


M 24 


244823040 


4.967107. . . 


29.84845. . . 


5p(4, F 7 ) 


276595200 


3.501127... 


14.83811. . . 


^ + (4, F 3 i) 


442828800 


3.829841. . . 


17.60003. . . 


fi~(4, F 3 i) 


443751360 


3.003133. . . 


11.02613. . . 




696729600 


4.194248. . . 


20.79438. . . 


McL 


898128000 


4.561453. . . 


27.45649. . . 


Sp(A, F 9 ) 


3443212800 


3.409108. . . 


14.04475. . . 


e 


4UoUoo 1 zUU 


o.4ooboU. . . 


i/l 011 in 
14.olliy. . . 


G 2 (F 5 ) 


5859000000 


3.855868. . . 


18.68766... 


^P(6,F 3 ) 


9170703360 


3.871692... 


18.90072... 


Co 3 


495766656000 


4.535119... 


25.99974. . . 


Co 2 


42305421312000 


3.865290. . . 


17.74829. . . 


fi(5,F 3 i) 


409387254681600 


3.277801... 


12.90986... 


Rub 


43252003274489856000 


5.668645. . . 


36.78701... 
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